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Abstract 

We study the feasibility of detecting the stau neutralino (fi-%?) coannihilation region at 
the LHC using tau (r) leptons. The signal is characterized by multiple low energy r leptons 
from X2 ~^ ~~ > TT Xi decays, where the fi and Xi mass difference (AM) is constrained to 
be 5-15 GeV by current experimental bounds including the bound on the amount of neutralino 
cold dark matter. Within the framework of minimal supergravity models, we show that if 
hadronically decaying r's can be identified with 50% efficiency for visible pr > 20 GeV the 
observation of such signals is possible in the final state of two r leptons plus large missing 
energy and two jets. With a gluino mass of 830 GeV the signal can be observed with as few 
as 3-10 fb _1 of data (depending on the size of AM). Using a mass measurement of the r pairs 
with 10 ftT 1 we can determine AM with a statistical uncertainty of 12% for AM = 10 GeV 
and an additional systematic uncertainty of 14% if the gluino mass has an uncertainty of 5%. 



1 Introduction 



The advent of supersymmetry (SUSY) has led to the construction of particle theory models 
linking a remarkably wide range of physical phenomena. While initially proposed on anesthetic 
grounds that nature should be symmetric between fermions and bosons, the fact that super- 
symmetry allows for the cancelation of the quadratic Higgs divergence allows the building of 
consistent models valid up to the grand unification (GUT) or Planck scales. The extension 
of supersymmetry to a local gauge theory, supergravity [TJ |2] , led to the development of GUT 
models giving a description of physics from Mqut down to the electroweak scale and incorporat- 
ing the successes of the Standard Model (SM) jHlSE]- Subsequently, LEP data experimentally 
confirmed the validity of the idea of SUSY grand unification. Further, an additional feature 
of SUSY is that models with i?-parity invariance automatically give rise to a candidate, the 
lightest neutralino (xi), for the astronomically observed cold dark matter (CDM) deeply link- 
ing particle physics with cosmology, and detailed theoretical calculations 6J confirm that GUT 
models can also achieve the experimentally observed amount of dark matter [7j in a natural 
way. A large number of experiments are now under way to try to detect SUSY dark matter in 
the Milky Way. Thus it is possible to construct models that encompass the full energy range of 
particle physics and simultaneously reach back into the early universe at times 10 -7 seconds 
after the Big Bang. 

If supersymmetry ideas are correct, then the Large Hadron Collider (LHC) is expected to 
produce the Xi particles. Similarly, the large one ton dark matter detectors, currently under 
development, will cover a large amount of the allowed SUSY parameter space and may detect 
directly the Milky Way dark matter particles. Theory predicts that these two should be the 
same, and the question arises as to how one might verify this. The direct approach to this 
problem would involve having the dark matter detectors measure the nuclear differential cross 
section for the incident dark matter particle scattering from the detector atomic nuclei, and 
compare these with the differential cross sections occurring for the neutralinos produced in the 
LHC. Such measurements of course may be many years in the future, and one would like to 
see if more immediate measurements might give strong indications for the equivalence of the 
astronomical and accelerator phenomena. To investigate this it is necessary to chose a SUSY 
model, and we consider in this paper mSUGRA [3J H] (which has been used in many other LHC 
calculations). 

The allowed mSUGRA parameter space, at present, has three distinct regions jB]: (i) the 
stau neutralino (^i-Xi) coannihilation region where Xi is the lightest SUSY particle (LSP), 
(ii) the Xi having a larger Higgsino component (focus point) and (iii) the scalar Higgs (A , H°) 
annihilation funnel (2M~o — M40 #0). These three regions have been selected out by the CDM 

Al 

constraint. (There stills exists a bulk region where none of these above properties is observed, 
but this region is now very small due to the existence of other experimental bounds.) The 
distinction between the above regions can not be observed in the dark matter direct detection 
experiments. It is therefore important to investigate whether the dark matter allowed regions 
can be observed at the LHC where the particles will be produced directly and their masses will 
be measured. The three dark matter allowed regions need very precise measurements at the 
colliders to confirm which is correct. 

In this paper we choose to work with the f\-x\ coannihilation region at the LHC. We 
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note that many SUGRA models possess a coannihilation region and if the muon magnetic 
moment anomaly, a^, maintains, it is the only allowed region for mSUGRA. Coannihilation is 
characterized by a mass difference (AM) between f\ and Xi of about 5-15 GeV in the allowed 
region. This narrow mass difference is necessary to allow the fi's to coannihilate in the early 
universe along with the xV s i n order to produce the current amount of dark matter density of 
the universe. Thus if this striking near degeneracy between f\ and Xi is observed at the LHC, 
it would be a strong indirect indication that the Xi was the astronomical dark matter particle. 
The coannihilation region has a large extension for the gaugino mass m^, up to 1-1.5 TeV, 
and can be explored at the LHC unless tan/3 is very large. 

We show here the feasibility of detecting the signal from the coannihilation region in the 
first few years of LHC running. The main difficulty, however, in probing this region is the small 
AM value. Staus from f\ — > TXi decays generates signals with very low energy tau (r) leptons 
and thus makes it difficult to discover this region at any collider due to the large size of the 
background events from the SM and SUSY processes. Signals with r leptons in the final states 
arising from the fi, x!ii an d Xi decays have been studied before in the context of the LHC for 
a large AM value (~50 GeV) |H]. A detailed Monte Carlo (MC) study has been performed 
to determine SUSY masses by fitting various mass distributions in SUSY particle decays [H], 
but at a pure generator level without any detector simulation such as effects by reconstruction, 
mis- identification, and the SM background. Those studies do not answer whether we can detect 
the coannihilation region at the LHC. The primary feature of this region is that the branching 
fraction of the staus and gauginos to r leptons is very large, and that the kinematics of small 
AM values produce these r's with very low energy. Both conspire to make any analysis difficult 
as there can be very large backgrounds for low energy r final states. However, they do suggest 
possibilities for measurements. Previously we have shown the procedure to measure small AM 
at an International Linear Collider (ILC) with an accuracy of about 10% for our benchmark 
points [TOl- But the ILC is still in the proposal stage. Thus it is important to see if the 
coannihilation signal can be observed at the LHC. 

We describe here an experimental prescription for detecting SUSY signal in the coannihi- 
lation region and measuring small AM at the LHC which provides a necessary condition to 
be in the coannihilation region. We first discuss the available mSUGRA parameter space in 
Sec. 2 and the sparticle masses in this region. In Sec. 3 we detail a search for r pairs from 
X° —* t?\ — > ttxi decays and select events in the two r leptons plus large jet(s) and large 
missing transverse energy (fir) final state. In particular, we motivate our selection criteria, de- 
scribe the tt mass and counting observables that would provide consistency checks and other 
evidence, and indicate the amount of luminosity needed to establish a 5a significance. In this 
section it will become clear that a primary consideration for this analysis to have sensitivity 
is for the LHC detectors to have efficiency for r's with > 20 GeV. In Sec. 4, we discuss a 
method of measuring AM and show the potential accuracy. We conclude in Sec. 5. 

2 mSUGRA Model in the Coannihilation Region 

The mSUGRA model depends upon four parameters in addition to those of the SM. They are 
mo, the universal scalar soft breaking parameter at M GU t; ^1/2, the universal gaugino mass 
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Figure 1: The narrow AM coannihilation band is plotted as a function of vnxii for tan/3 = 40 
with A = and \i > 0. The left end of the band is due to the b — > 37 branching ratio bound 
and the right end by < 11 x 10~ 10 . 

at Mgut; A , the universal cubic soft breaking mass at Mqut; and tan/? = (Hi) / (H 2 ) at 
the electroweak scale, where Hi (H 2 ) gives rise to up-type (down-type) quark masses. The 
model parameters are already significantly constrained by different experimental results. Most 
important for limiting the parameter space are: (i) the light Higgs mass bound of M^o > 
114 GeV from LEP [TTj . (ii) the b —> 57 branching ratio bound of 1.8 x 10~ 4 < B(B —> 
Xgj) < 4.5 x 10 -4 (we assume here a relatively broad range, since there are theoretical errors 
in extracting the branching ratio from the data) ^2], (hi) the 2a bound on the dark matter 
relic density: 0.095 < f2 CDM /i 2 < 0.129 [7j, (iv) the bound on the lightest chargino mass of 
M~± > 104 GeV from LEP and (v) the muon magnetic moment anomaly a M , where one 

gets a 2.7a deviation from the SM from the experimental result El El- Assuming the 
future data confirms the a M anomaly, the combined effects of g M — 2 and M~± > 104 GeV then 
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only allows /i > 0. 

Figure ^ shows the range of allowed AM values in the coannihilation region as a function 
of mi/2 for tan/3 = 40. We see that AM is narrowly constrained and varies from 5-15 GeV. 
Because of the small AM value, r's from fi — > txi decays are expected to have low energy 
providing the characteristic feature of the coannihilation region. From here on, we assume 
m i/2 — 360 GeV and tan/3 = 40 as a reference point and use it throughout the text unless 
otherwise noted. SUSY masses for five example those points are listed in Table fusing ISA JET 
7.63 PH. 
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Table 1: Masses (in GeV) of SUSY particles in five representative scenarios for mi/j = 360 GeV, 
tan/? = 40, [i > 0, and A = 0. These points satisfy all the existing experimental bounds on 
mSUGRA. M^. ax is the end point of the di-tau mass distribution from the x® decays. The 
production cross section is nearly same for the five points and is 8.3 pb for ttlq = 215 GeV. 
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3 Detecting a SUSY Signal in the Coannihilation Region 

Unlike at the ILC ^Uj, it is very difficult to access direct r^rf production at the LHC because 
of its small production cross section relative to the huge QCD backgrounds. Therefore, we will 
show that using copious X2 an d decay via X2 ~^ T ^1 "~ > T T Xi a significant number of di-tau 
events can be observed with an effective observable. The end point of the mass distribution is 
calculated to be 



M T m T ax = M~o, 
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which corresponds to the case when the two r's are back-to-back in the X2 Tes ^ frame. To 
reconstruct the mass as fully as possible we use hadronically decaying r's (rt's). Despite the 
fact that the r's lose energy to neutrinos, the di-tau mass distribution (M y f) still provides a 
peak position (M^ ak ) that can be directly measured. The peak position is below M^. ax and 
M™ is effectively smeared. We next describe how to use these experimental features to detect 
the signal. 

The primary SUSY production processes at the LHC are pp — > qg,qq,gg- In each case the 
decays proceed via q — > q'xf or q\\ (or q R — > g ~^ qq'xf or qqx\'i an( ^- 9 ~^ tti or ^1 ari d 
their charge conjugate states, generally producing high E T jets and gaugino pairs. We are most 
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interested in events from X1X2, X1X21 or X2X2 P arrs ; where the Xi m the first case is directly 
from the <jr decay. The branching ratio of X2 ~* T ^1 i s about 97% for our parameter space 
and is dominant even for large m^ 2 in the entire coannihilation region; the same is true for the 
Xi — > v f\ decay mode. (It should be noted that both e# and \xr are lighter than x\ by about 
10 GeV . However, the branching ratio for x\ — > eeR or fifiR is much less than 1%.) Since the 
stau decays via fj — > TXi, we expect inclusive X2 events to include at least two r leptons plus 
large E T jet(s) and large fir (from the xT)- 

We consider two experimental scenarios. The first uses the fir + >2 jet final state to reduce 
backgrounds and searches for the 2r's that arises from the decays of yQ. In each candidate 
event, all di-tau pairs can then be searched for a mass peak as evidence of the X2 decay chain. 
The second option studies gaugino pairs (X2X2; X1X2) an d requires 3r's to reduce backgrounds 
and only 1 jet + fir in the final state to regain acceptance. Both final states will be triggered by 
requiring large E T jet(s) and large fir and such a trigger will be available at both the ATLAS 
and the CMS experiments. We focus on the final event selection and report here on an analysis 
of the fir + >2 jet + > 2r events using the ISAJET event generator and the ATLAS detector 
simulation ATLFAST 18 j. The results of the 3r analysis will be considered elsewhere |19j . 

To establish a signal in the coannihilation region we require (at least) (1) a 5a excess in the 
counting of + 2 jet + 2r events with M™ below M^. ax and (2) the mass peak to be consistent 
with the expectations from X2 decays. The question then becomes whether such measurements 
are possible at the LHC. We begin by addressing the primary experimental issues: (a) the p-r 
spectrum of r's from fi — > tx\ which is expected to be soft and (b) the ability to select the 
correct di-tau pairs. 

We first examine the visible pr (p™) distribution of Th's from f\ — > TXi in X2 decays with 
AM = 5, 10, and 20 GeV using ISAJET. As seen in Fig. El even with a small mass difference, the 
r is boosted in the cascade decay of the heavy squark and gluino making it potentially viable. 
One can already begin to see the importance of reconstructing the r's with p™ > 20 GeV. From 
here on, we assume that both the ATLAS and CMS detectors can reconstruct and identify Th's 
with p™ > 20 GeV. Assuming the taus are identifiable, we require two jets with Et > 100 GeV 
and fir > 180 GeV. These cuts should satisfy the fir + jets trigger at ATLAS and CMS with 
high enough efficiency that we further neglect trigger efficiency and bias effects. 

Events are further selected by requiring at least two identified t^s with at least one r 
lepton with p^ s > 40 GeV (from x° ~~ > r ^i) an d one with p™ > 20 GeV (from f% — > TXi)- At 
this point, the dominant background is expected to be from ti pair production, where each t 
decays as t —>■ bru and produces a final state of fir + 26 + It. A correlation plot between fir and 
£ T jetl +£ T jet2 is shown in Fig.Hfor ti and SUSY events. We choose £ T jetl + £ T jet2 + fir > 600 GeV 
to further reduce the background. 

For each candidate M^ s is calculated for every pair of r's in the event and categorized as 
opposite sign (OS) or like sign (LS). The mass distribution for LS pairs is subtracted from 
the distribution for OS pairs to extract x\ decays on a statistical basis. The result is shown in 
Fig. 13] where we have assumed the identification (ID) efficiency (e) to be 100% and a probability 
that a jet is misidentified as r h (/j_> T ) to be 0%. We see that the non-%2 OS pairs are nicely 
canceled with the wrong LS combination pairs and that the OS— LS distribution is well fit to 
a Gaussian. We note that while the expected maximum di-tau mass, M^. ax = 78.7 GeV, is 
consistent with Fig. 0] and easily determined, we believe a full detector simulation is required 



5 



> 

CD 

u 

CM 
W 

c 

3 
O 

o 



-Q 

E 



i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i | i i i i I 

(from %° decays) 

— AM = 5GeV 

— AM = 10 GeV 
AM = 20 GeV 




10 20 30 40 50 60 70 80 90 100 

Visible P t (t) (GeV) 

Figure 2: The visible distributions for hadronically decaying r leptons from f\ — > 7~x5 where 
the fi's are required to be decays of x° — * r ^i- The curves show the results for AM = 5, 10, 
and 20 GeV. 
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Figure 3: Scatter plots of fa vs E^ etl + E^ et2 for tt (left) and SUSY (right) events after 
requiring 2r's passing all ID and kinematic cuts in the final state. The SUSY events are for 
our reference point at M~ g = 831 GeV and AM = 10.6 GeV. We require £ T jetl > 100 GeV, 

£ T jet2 > 100 GeV, fa > 180 GeV, and £ T jetl + £ T jet2 + fa > 600 GeV to eliminate the top 
background. 
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Ditau Invariant Mass (GeV) 

Figure 4: The visible di-tau invariant mass (M^ s ) distributions for all di-tau combinations in 
inclusive fir + 2j + 2r events in the coannihilation region (AM = 10.6 GeV; see Table ^) 
with e = 100% and fj-+ T = 0%. The r's are selected with p™ > 20 GeV but requiring at least 
one r lepton to have p™ > 40 GeV. The endpoint is consistent with the theoretical value at 
M™ ax = 78.7 GeV. The dashed and solid open histograms are for OS and LS pairs, respectively. 
The gray (or yellow) histogram is for OS— LS pairs and is well fitted to a Gaussian (solid line). 
The small excess between 80 and 150 GeV is mainly due to x° (456 GeV) — > rfi, ZyQ, W ± xf 
decays. 

to show this in practice. We note for completeness that there is a small second excess between 
80 and 150 GeV, which is mainly due to \\ (456 GeV) — > rfi, Z\% W ± xf decays. 

We next assume e = 50% and /j_> T = 1%, based on the CDF measurement [20], and show 
that these results are not significantly affected by these experimental conditions except to 
quadruple the needed luminosity to establish a signal. Figure El shows the M™ distributions 
for OS, LS, and OS— LS pairs for a hypothetical 10 fb _1 sample to allow visual confirmation of 
the peak taking into account statistical fluctuations. Since the number of it events that survive 
our selection cuts is expected to be a few OS— LS counts per 10 ftT 1 we do not consider them 
further. The excess between 80 and 150 GeV seen in Fig0] becomes statistically insignificant. 
It should be noted that a set of similar kinematical cuts is also used in Ref. |H] and the SM 
background was also found to be small. To further illustrate the importance of reconstructing 
r's with £>t 1s ~ 20 GeV, the RHS of Fig. E] shows that the mass peak disappears for the same 
sample, if we require all r's to be above 40 GeV. 

For any coannihilation signal to be established a necessary (but not sufficient) condition is 
that an excess of signal events above background should be readily observed. Figure El shows 
the expected yield of OS— LS counts below the end point per 10 fb _1 as a function of AM. 
For comparison, we show the results for e = 50% and fj-> T = 1%, along with the ideal r 
ID condition of e = 100% and /j-> T = 0% but scaled down by a factor 4 due to 1/e 2 . It is 
remarkable to see that there is no significant difference between two cases, suggesting that the 
counting is very insensitive to the fake rate. This suggests that the counting nearly perfectly 
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Ditau Invariant Mass (GeV) Ditau | nvariant Mass (GeV) 

Figure 5: [left] Same as in Fig. HI but using a randomly-selected 10 fkT 1 SUSY sample that is 
analyzed with e = 50% and fj^ T — 1%- M^. ax = 78.7 GeV can still be inferred and a clear peak 
is visible. Note that the it contribution in M™ < 100 GeV is estimated to be a few OS— LS 
counts and has not been plotted, [right] The same sample but requiring all r candidates to 
have p™ > 40 GeV. The characteristic di-tau mass distribution has disappeared. 



reflects the counting of x^'s if one takes into account the r ID efficiency This also indicates that 
in the signal region there is only a small contamination from the wrong combination events. 
FigureElalso shows the numbers of OS— LS counts if the gluino mass is varied by approximately 
±5% (and thus also the x[> an d Xi masses), but maintaining the same AM values. Heavier 
(lighter) gluino mass is reflected as a smaller (larger) production cross section, thus with a 
smaller (larger) yield. 

We next determine the luminosity necessary to observe OS— LS counts with a 5a signifi- 
cance. We define a = (N os - N LS )/y/N os + N LS , again require M™ < M™ ax , and use the 
results in Fig. El to parameterize the number of OS— LS counts as a function of AM j2TJ. The 
luminosity required for as a 5a excess as a function AM is shown in Fig We conclude that 
the characteristic coannhilation signal (AM = 5-15 GeV) can be detected at the LHC with as 
few as 3 fb" 1 if AM is large (~15 GeV) and with 10 fb" 1 if AM is small (assuming the SUSY 
GUT models are correct and the scale is Mg = 831 GeV). 

The second criteria we consider to help establish the coannihilation region signal is the 
observation of M^r ak at the correct location. Within mSUGRA models, since the g, \\ an d X® 
masses are related, M^ ax changes with M= for a given M~o. This also means M>^ ak changes 

1 Al 

with AM. Figure |H1 shows Mj^ ak as a function of AM; any experimental observation would 
need to be consistent with this small range of 39-52 GeV for the case of Mg = 831 GeV. 

We next comment on two effects that could that could change the expected M^: ak value. 
The first is that since there is a direct relationship between M~o and Mx, from Eq. d the 

peak position is expected to be a function of Mg for a given AM. If Mg is changed from 

our nominal value but the AM value is maintained, we find <5^p ea k /M^ ak ~ £ j^max /M™ ax = 

0.718 5m~/Mq. This is directly observed in Fig. Hand effectively expands the observation range 
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to 38-54 GeV. We also study the potential impact of varying r ID conditions. On the same 
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Figure 6: The mean number of OS— LS counts below the end point per 10 fb" 1 of data as a 
function of AM. The solid circles are for e = 50% and /j_> T = 1%. The solid line is the case 
with e = 100% and /j-> T = 0%, but with the number of OS— LS counts scaled down by a factor 
of 4 due to 1/e 2 . The uncertainty of each circle is about its own size based on 200 fb" 1 SUSY 
MC samples. The dashed and dash-dotted lines are when the gluino masses are varied by about 
5%. 




5 10 

AM [GeV] 

Figure 7: The luminosity necessary to establish a 5cr significance in the number of OS— LS 
counts for the coannihilation region. The band reflects a variation due to the gluino mass by 
±5% from a nominal gluino mass of 831 GeV in our reference model (see Table [Q). 
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Figure 8: The peak of the M™ distribution as a function of AM. The solid circles are for 
e = 50% and fj^ T = 1%, compared to the solid line for the case with e = 100% and fj^ T = 0%. 
The statistical uncertainty on each circle is about the size of the circle from 200 fb _1 SUSY MC 
samples. The dashed and dash-dotted lines are when the g masses are varied by about 5% and 
can be understood as coming from the change in the X2 m ass for a fixed AM, since mSUGRA 
provides a fixed relationship between Ms and M~o (See Eq.l). Note that the variation due to 

" X2 

the fake rate is negligible. 



figure we show (a) e = 100% and fj^ T = 0% and (b) e = 50% and /j-_» T = 1% and see that 
there is no systematic peak shift due to the jet — > r fake rate. This can be explained by the 
fact that the jet misidentified as a r has no correlation with the r leptons from the x\ decay; 
both OS and LS di-tau mass distributions are affected equally in its shape, and thus will cancel 
in the difference to good approximation. 

As an aside we note that there is a non-coannihilation region where a low energy r lepton 
is expected when the fi mass is close to the x\ mass instead of the Xi- 111 this "inverted mass" 
scenario, AM inv = M~o — M~ = 5-15 GeV. SUSY masses in two particular cases are listed 

A 2 1 

in Table |2l compared to our reference point (AM = 10.6 GeV and M^. ax = 78.7 GeV). The 
di-tau mass distribution for the events would look similar to those in the coannihilation region, 
but we expect a difference in the invariant mass of the jet-Th-Th system in the q — > qx\ decay 
where the squarks are heavier by about 30 GeV compared to those of the reference point. A 
full study is beyond a scope of this study, and we will address this in a future paper. We note, 
however, that the above two cases do not satisfy the relic density constraint and hence are 
inconsistent with cosmology. Also, while the small AM of the coannihilation region is natural 
for mSUGRA to be consistent with WMAP data, the inverted mass scenario further requires 
an artificial fine tuning to make M~o — M~ small. 

A2 1 

Having shown that we can establish a 5<r excess and that the peak and end point of di-tau 
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Table 2: Masses (in GeV) of SUSY particles in two "inverted mass" scenarios for mj/2 = 
360 GeV, tan/3 = 40, /i > 0, and A = 0. 



m 323.7 315.2 

~& L 798 795 

*i 598 560 

X° 2 266.0 266.0 

fi 255.4 247.8 

X? 144.2 144.2 

= M~o-M fi ) 10.6 18.2 

M™ ax 61.4 78.7 



mass distribution are consistent with expectations and only occupy a small region of possible 
values, we next discuss the prospects for measuring AM with 10 fb _1 of data using the peak 
position of the di-tau mass distribution. 

4 Measuring AM Using 10 fb 1 

The pr+ 2j + 2r data set provides an opportunity to make a AM measurement. Since 
varies with AM, a measurement of M^r ak combined with a measurement of Ms performed 
elsewhere allows for a measurement of AM. It should be noted that from Fig. [HI the peak 
position varies with M%. This arises from the mSUGRA relations M% ~ 3.1 M~o and the fact 

y y A2 

that the gaugino masses enter in the endpoint relation of Eq. ^ We note for completeness that 
a similar measurement could be done with the results of Fig. |3 This will be done is a separate 
paper [T9"j . 

The statistical uncertainty of the AM measurement is dominated by the precision of the 
Mpeak measurement. The statistical uncertainty in the measurement of MJ^ ak is typically in the 
1 GeV range for 10 fb" 1 . Numerically it is given by um/ y/N, where (Ju is the r.m.s. of the mass 
distribution (see Fig. HJ) and is of the order of 10 GeV, and N is the number of OS— LS counts 
in the peak region for a given luminosity (see Fig. |BJ). For a more detailed parameterizations 
see Ref. j2TJ. The systematic uncertainty due to the uncertainty in the assumed value of Mg 

which changes the M^^ vs AM relationship (see Fig. |HJ) is also of the order of 1 GeV for 
a 5% uncertainty in the gluino mass. The uncertainty due to fj^ T does not affect the peak 
measurement. 

For the cases of AM = 6 and 10 GeV, with 10 fb~ of data, we estimate that we can 
determine the AM values to be 



AM 




V9 GeV 
±1 GeV 
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Figure 9: The uncertainty in the AM measurement as a function of luminosity. The solid 
(dashed) curve indicates where this analysis has more than (less than) 5a significance in OS— LS 
counts and shows the expected statistical uncertainty. We also show the systematic uncertainty 
(horizontal dotted line) in AM due to a gluino mass uncertainty of ±5%. 



where the first uncertainty is due to the statistics and the second due to a gluino mass variation 
of 5% from the nominal value of 831 GeV. If one assumes a gluino mass uncertainty of 10%, the 
second uncertainty is almost doubled to tf/? and j|; 3 for the 6 and 10 GeV cases, respectively. 
Figure M illustrates both uncertainties for the AM = 10 GeV case as a function of luminosity. 



5 Conclusion 

In mSUGRA models, the fi-Xi coannihilation occurs in a large region of parameter space 
allowed by the relic density constraint. The characteristic feature of this coannihilation region 
is a small mass difference (AM) between the f\ and the Xi- The small mass difference produces 
final states containing low energy r's from x\ — > rf\ — > ttXi- We have demonstrated that if 
LHC experiments reconstruct /identify r's with > 20 GeV with an efficiency in the 50% range, 
we could establish the signal for this coannihilation region using a sample of fir+ 2j + 2r events. 
Counting the number of OS— LS di-tau pairs in conjunction with the observation of a mass peak 
below the expected end points would establish the signal for the case of Mg ~ 830 GeV with 

as few as 3-10 fbT 1 of data. We expect to require a larger luminosity as Mg increases. We 
note that although we have assumed that the jet — > r fake probability was 1%, our results 
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appear to be insensitive to the fake rate. Using the di-tau mass measurement and assuming 
the gluino mass is already measured, we have further shown that a measurement of AM is 
possible even with 10 fb _1 of data. For AM = 10 GeV, the statistical uncertainty is 12%, with 
a 14% additional uncertainty due to an assumed 5% uncertainty in the gluino mass. Again, 
the results are insensitive to the fake rate. 

We note that the signal of $t+ ljet + 3r also occurs at a reduced rate from gaugino pairs, 
but with lower backgrounds. However, it can provide a complementary signal that could help 
solidify the discovery of the coannihilation signal at the LHC. 

The coannihilation region is also present in most non-universal SUGRA models, and an 
analysis similar to the one here for mSUGRA could be performed for these. The main require- 
ments for an observable signal would be that the x° decays predominantly to r's, and that the 
decay of the parent gluinos and squarks to X2 * s n °t suppressed. 
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